Time evolution of an anharmonic oscillator interacting with a squeezed bath 

G. M. D'ArianoQ 

Dipartimento di Fisica "Alessandro Volta", Universita di Pavia, Via A, Bassi 6, 1-27100 Pavia, Italy 

M. FortunatdO 
Dipartimento di Fisica, Universita di Roma "La Sapienza" , P.le A. Mora 2, 1-00185 Roma, Raly 






o 



and 

P. Tombesi 
Dipartimento di Matematica e Fisica, Universita di Camerino, Via Madonna delle Carceri, 1-62032 Camerino, Raly 

(Received February 1, 2008) 



The evolution of a single mode of the electromagnetic field interacting with a squeezed bath in a Kerr medium is considered. 
The solution of the corresponding master equation is given numerically. It is argued that the creation of a superposition state 
(Schrodinger's cat) is better achieved in presence of a squeezed reservoir than of a thermal one. 
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I. INTRODUCTION 

The possibility of detecting quantum interference ef- 
fects between macroscopically distinguishable states — 
eigenstates of a variable corresponding to macroscopi- 
cally different eigenvalues — has received much attention 
in the last few years [0-|O|. Such interest is motivated 
by one of the fundamental problems in quantum mechan- 
ics: to determine whether quantum-mechanical features 
can be observed in macroscopic objects ||ll[. Unfortu- 
nately, the detection of a macroscopic superposition is 
very difficult, due to the unavoidable coupling with the 
environment and by the consequent dissipation |l3] . The 
severity with which dissipation destroys the quantum co- 
herence at macroscopic level has been discussed by sev- 
eral authors ||l^-[l^. 

Since the original proposal by Yurke and Stolcr ffl , who 
have shown that a coherent state propagating through 
a Kerr medium evolves, under suitable conditions, into 
a quantum superposition of two coherent states which 
are 180° out of phase with each other, many attempts 
have been made in order to preserve this superposition 
with detectable effects as, for example, the interference 
fringes at the output of a homodyne detector. Mecozzi 
and Tombesi |3|-|g] have considered a dissipation model 
corresponding to a beam splitter with a squeezed vac- 
uum injected into the unused port. They have shown 
that the interference pattern may be preserved if the 
input light is squeezed in a suitable quadrature. As a 
model for phase-sensitive measurements, Kennedy and 
Walls [|| suggested the use of a squeezed bath in place of 
the thermal one, showing a substantial improvement of 
the quadrature-phase sensitivity. On these lines, Buzek, 
Kim and Gantsog |1^ have recently studied the phase 



properties of superpositions of coherent states evolving 
in a squeezed (phase-sensitive) amplifier, showing that 
differently from the phase-insensitive case [|9| , the phase 
distribution of the Schrodinger-cat input state can be 
preserved for long times. 

The above models either suffer a parametrically im- 
posed dynamics for the loss (beam splitter models 
H'ItIIDj 01' do not consider the nonlinear coupling gen- 
erating the superposition state, which is just treated as 
an initial condition p,pJ3|. For this reason we consider 
a model in which both the anharmonic Hamiltonian and 
the interaction with a squeezed bath are taken into ac- 
count, thus allowing a test of the effect of squeezed fluctu- 
ations on the generation of a superposition state starting 
from a single coherent one. 

This article is organized as follows: in Sec. O we in- 
troduce the model and discuss the corresponding master 
equation. In Sec. Ill we present the numerical integration 
of the master equation and the results of the numerical 
analysis, showing the efficiency of a squeezed bath in the 
creation of a "quasi-superposition state" . In Sec. IV we 
conclude and summarize the results. 



II. THE MODEL AND THE MASTER EQUATION 

We consider a mode of the electromagnetic field at fre- 
quency Lu in a Kerr medium coupled to a squeezed bath, 
namely a reservoir of oscillators whose fluctuations are 
squeezed H. The total Hamiltonian is given by 



H — Hs + Hi + Hb , 



(1) 



where iJa is the free Kerr Hamiltonian of the field mode 



Hs = ujia^^a) + n{a^a)^ 



(2) dQ 



H-Q is the free Hamiltonian of the bath, and Hi is the 
osciUator-reservoir interaction Hamihonian, which in the 
rotating wave approximation has the form 



H, = a)t + aft 



(3) 



In Eqs. (||) and (||) a^ and a are the boson creation and 
annihilation operators of the mode, while F^ and F are 
bath operators. The bath is squeezed and Markovian, 
namely the correlations functions of the operators F^ and 
F are given by Q 

(ft(i)f(t'))=27iV<5(t-t'), 
(f(t)ft(t'))-27(iV+l)5(t-0, 
(f (t)f(t')) = 27Af e-2^'^*(5(t - t') , 
(ft (i)f t (t')) = 27M*e2*'^*,5(t ~ t') . (4) 

In Eqs. (U) 7 is the damping constant (determined by the 
coupling between the oscillator and the bath) , A^ is a real 
parameter, which reduces to the mean number of thermal 
photons when the bath is not squeezed, and M is the 
squeezing complex parameter {M = |Af je*''^) satisfying 
the relation 



\M\ < y/N{N + 1) . 



(5) 



The parameters N and M measure the strength of the 
correlations of the bath degrees of freedom: for |Af p = 
N{N+1) the squeezing is maximum, whereas for |il/| = 
the reservoir reduces to the customary thermal one. 

From Eqs. (|ll)-(|j) one can derive the following master 
equation for the reduced density matrix of the field mode 
in the interaction picture |20 



-^ = -in[{a^af,p] +7(iV + l)(2apat - a'' ap - pa'' a) 



+ ^N{2a* pa — aa* p — paa*) 

- 7M(2at/5at - a'^a'^p - pa'' a'') 

— jM* (2apa — aap — paa) . 



(6) 



Eq. (ph is numerically solved in the next section. The 
results are shown in terms of the Q-function 



Q{a,a\t)^{a\p{t)\a) , 



(7) 



which is the (positive definite) probability density for the 
antinormally ordered moments of the annihilation and 
creation operators, and in terms of the Wigner function, 
which is defined by 



W{a,a*,t)= I ^e 



-Aa*+A*Q 



Tr 



{p(i)e^"'-^*«} . (8) 



Using standard methods |1^, it is possible to con- 
vert the master equation (0) into the Fokker-Planck type 
equation for the Q-function 



dt 



= -in 



a*(l + 2|ang 



+ 7 



d 



cPQ_ 
d 



a(l + 2|aP)g 






a 



da 



d^Q 



^^,«* g + 27(A^+l)-^^^^, 



■7M 



.d^Q 



(9) 



Eq. (|) has been solved for f7 = [| and for M = 0^. 
For ri and M both nonvanishing, however, an analytic 
solution is not available, whereas numerical integration 
has to be carried out carefully, because of computation 
instabilities pi[] . 

In absence of dissipation an initial coherent state 



ao 



exp 



\ao\ 






n=0 



(10) 



evolves towards the following superposition of coherent 
states at t = 7r/251 



3=(e-/Vo) 



V2 



t/4| 



ao, 



(11) 



The Q-function of the t = coherent state is the Gaus- 



Q(a, a*,0) = exp(-|a- aop) , 



(12) 



whereas for t — 7r/2J7 it approximatively corresponds to 
two Gaussian peaks in the complex plane centered at ao 
and — ao- The two peaks are "macroscopically distin- 
guishable" for |ao| ^ 1. The evolution is periodic with 
period T = 27r/r2. In presence of dissipation the state 
at times t > is no longer pure, and a superposition 
of states corresponds to a "quasi-superposition" mixed 
state, with interfering features still surviving. We will 
look for signatures of such states from peaks in the com- 
plex plane approximately separated by a distance 2|ao|, 
reminiscent of the state (pTl). 



III. NUMERICAL RESULTS 

We have studied the time evolution of the density ma- 
trix by integrating numerically the master equation m) 
for truncated Hilbert space dimension d = 128 (typically 
a power of 2, in order to take advantage of fast Fourier 
transform algorithm for reconstructing the Q-function 
and the Wigner function). The integration time-step has 
to be carefully tuned as a function of Q and ao (here 
typically jAt ~ 10~^ for a standard fourth-order Runge- 
Kutta routine). Numerical accuracy is checked through 
normalization of /5, Q and W , positivity of Q, and reality 
of the diagonal elements of p. As a test, the results from 
analytical solutions for 7 = ||3,0] and for fJ = [pi have 
been recovered up to the seventh digit. 



In Fig. H the Q-function from the maste r equation 
(|) for maximaUy squeezed bath [|M| = y/N{N + 1)] 
is given for il/j — 10, N — 30, and squeezing phase 
■ip = 0. At 7i = 0.0576 [well before fit = 7r/2, namely be- 
fore the creation time of the state (O) in the undamped 
case] two peaks are clearly visible in the structure of the 
Q-function. This result can be compared with that of 
Daniel and Milburn |1^ at the same time, but in absence 
of squeezing. The two peaks are almost symmetrical with 
respect to the origin of the phase space. Of course, due 
to dissipation, they are not strictly coherent, but they 
are still visible, even for subsequent times. 

In Fig. the effect of a purely thermal bath on the 
generation of the superposition state is shown for com- 
parison with Fig. |l| for the same values of the parameters. 
In this case there are no peaks in the Q-function, which 
now exhibits a classical behaviour [|l^: this is the usual 
effect of dissipation |1J] . 

In Fig. H we show for completeness the time evolution 
of the Wigner function at the same times as in Fig. ^: 
the two peaks still show up, but now some interference 
features arise between the two component states, as ex- 
pected. We have also considered the case in which the 
phase of the complex squeezing parameter M is chosen 
as 4' — T^i namely the squeezing direction of the phase 
sensitive bath is rotated of 90° with respect to the case 
■0 = |lG,p2|: the creation of two peaks is still visible 
and their distance is approximately the same as in Fig. |l|, 
but the line connecting them is now rotated of about 90°, 
in agreement with the phase of Af . 

The effect is further confirmed by the behaviour of the 
quadrature distributions P{x^), defined by 



P{x^) = {x4,\p\x^) 



with 



X, = -^(e 



itp 



V) 



(13) 



(14) 



The probability P{x^) is plotted in Fig. | at 7* = 0.0576 
for = 7r/4 and = 37r/4. For = 37r/4 — along 
the direction which joins the two peaks of the Wigner 
function — the probability P{xcj,) exhibits two peaks re- 
sulting from marginal integration in the complex plane. 
On the other hand, along the orthogonal direction at (p — 
7r/4, a pattern reminiscent of interference in phase space 
^3| between the two originally coherent components is 
visible. We have also computed the degree of mixing 
S = 1 — Tr(p^): its behaviour shows a rapid increase with 
time, achieving the asymptotic value 5* = 0.96. Thus, 
the above coherence features still survive notwithstand- 
ing the large mixing. Notice that, however, this happens 
only for the two-peak state, whereas there is no reminis- 
cence of the symmetrical superposition of three coherent 
states which arises at t = 7r/3f2 in absence of dissipation. 
As regards the other parameters, the overall scenario is 
not affected by changing the phase of the initial coherent 



field ao: whereas increasing the modulus of ao deterio- 
rates the formation of the two peaks. On the other hand, 
slightly increasing f2 or N (still for a maximally squeezed 
bath) , does not appreciably improve the visibility of the 
peaks. 



IV. SUMMARY AND CONCLUSIONS 

We have considered the anharmonic oscillator inter- 
acting with a squeezed bath. The dissipative Kerr ef- 
fect models a coherent state propagating into a nonlin- 
ear medium, as, for example, into an optical fiber. The 
squeezing has been considered in order to have the pro- 
duction of coherent effects in presence of dissipation. Dif- 
ferently from the previous works |6| and ||l^, our model 
treats all effects (nonlinearity, dissipation and squeezing) 
contextually We have shown that squeezing the fluc- 
tuations of the bath improves the generation of "quasi- 
superposition" states, which are completely forbidden in 
presence of a heat bath. Our results are in agreement 
with those of Kennedy and Walls and of Buzek, Kim 
and Gantsog ||l^ and complete them: a squeezed bath 
not only is able to preserve a macroscopic quantum super- 
position, but also allows its generation in presence of non- 
linear interactions. The physical meaning of a squeezed 
bath is not completely understood in the present con- 
text: a suitable feedback mechanism could be envisaged 
that supports an ad hoc phase sensitive interaction, which 
amplifies or attenuates fluctuations, depending on the 
quadrature p^ . 

In conclusion, some comments regarding the direction 
of squeezing in the bath are in order. Despite the squeez- 
ing direction does not affects the generation of quasi- 
superposition states, the survival time of coherence is 
naturally increased only for squeezing in the direction 
orthogonal to the line that joins the component states. 
On the other hand, the Kerr effect rotates the Q-function 
in the complex plane, thus making squeezing less efficient 
in the overall evolution. This suggests that an ideal bath 
should be squeezed isotropically in the plane, depending 
on the phase of the state itself. It is possible to write a 
master equation for such an isotropically squeezed bath: 
numerical results on these lines will be the object of a 



forthcoming paper |25|. 
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FIG. 1. Contour plots of the time evol ved Q-func tion ac- 
cording to Eq. (|). Here N = 30, |M| = ^iV(iV-Hl), V = 0, 
ao = 3, and fl/'y — 10. 



FIG. 2. The same as in Fig. §, but for M = 0. 



FIG. 3. The Wigner function corresponding to plots in 
Fig. I 



FIG. 4. Plot of the marginal distribution P{x^) at 
■yt = 0.0576 for 4> = 37r/4 and = 7r/4. 
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Fig.l G. M. D'Ariano, M. Fortunate, and P. Tombesi 
Time evolution of an anharmonic oscillator ... 
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Fig.S G. M. D'Ariano, M. Fortunate, and P. Tombesi 
Time evolution of an anharmonic oscillator ... 
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Fig. 3 G. M. D'Ariano, M. Fortunate, and P. Tombesi 
Time evolution of an anharmonic oscillator ... 
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Fig. 4 G. M. D'Ariano, M. Fortunate, and P. Tombesi 
Time evolution of an anharmonic oscillator ... 



